Phase Diagram of a Two-Dimensional Neutral Classical 
Coulomb Gas from a Non-perturbative sine-Gordon Expansion 
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Within the sine-Gordon formalism of a two-dimensional neutral classical Coulomb gas, a conver- 
gent expansion with non-perturbative nature is performed to calculate the thermodynamic potential 
and construct the phase diagram. It is shown that truncation at the first order yields the Gaussian 
approximation. The second- and third-order corrections are analyzed for the case of small fugac- 
ity and are shown that they substantially improve the Gaussian-approximation phase diagram. In 
particular, these corrections introduce a new conducting phase and make the insulator-conductor 
coexistence phase end at the conducting phase. The latter result is in agreement with the prediction 
. made by generalized renormalization-group calculations. 

^ I PACS numbers : 05.70.Fh; 64.60-i; 64.60.Cn; 68.18.Jk . 

Two-dimentional neutral classical Coulomb gas (2DNCCG) can be mapped into a two-dimensional XY model [1|, 
and provides a prototype model for two-dimensional systems with which vortices are important thermal excitations [2|. 
These systems including supcrfluid films, superconducting films, Josephson-junction arrays, two-dimensional melting, 
surface roughening, liquid crystals and double-layer quantum Hall systems can undergo Berezinskii-Kosterlitz- 
C^l Thouless (BKT) phase transitions Q. Hence, the 2DNCCG has wide applications in studying these systems and, thus, 
O^l •& thorough understanding of its phase structure is of extreme importance. A number of analytical and numerical 
investigations were made to construct the phase diagram last three decades and they have greatly increased our 
understanding on the 2DNCCG. However, the picture is still far from complete and there remain several issues to be 
'settled as described below. 

Both the analytical and simulative investigations have revealed some universal features of the 2DNCCG. The Gas 
suffers a BKT phase transition when the fugacity or particle density is very small or low; that is, for a very samll 
fugacity, the gas remains as an insulator with a dipole phase when temperature is low, whereas, when temperature 
increases, the dipole charge becomes unbound, and, thus, the dipole phase becomes a conducting phase. This kind 
' <w> of phase transition is continuous. Moreover, increasing the fugacity from a small yet finite value, one encounters a 
first-order phase transition to an insulator-conductor coexistence phase. 

Nevertheless, the concrete phase diagrams for the 2DNCCG constructed by various methods have not reached a 
universal agreement. For example, the Kosterlitz's renormalization group (RG) equations at the lowest j| (1974) 
. £h an< ^ higher orders [|5| of an expansion in fugacity showed existence of just the BKT phase transion only, while 
K> Minnhagen's generalized RG equations [gj, Monte Carlo simulations generalized Debye-Huckel-Bjerrum theory 
j_j -|| , and Gaussian approximation [[l0| all led to both the BKT and the first-order phase transions. Monte Carlo 
Cd simulations gave even much richer phase structures M. For the connection between the BKT and the first-order 
phase transition curves, the Gaussian approximation ]10| showed that, in the fugacity- temper at ure plane, the first- 
order phase transition starts from the end point of the BKT phase transition curve ]10|| , and Debye-Hiickel-Bjerrum 
theory j^] produced the similar tricritical point in the density-temperature plane, whereas Minnhagen's investigation 
|(| and Monte Carlo simulations indicated that the first-order phase transition curve passes the end point of the 
BKT phase transition curve and ends at a higher temperature and in the conducting phase. Another point of dispute 
is the starting point of the BKT phase transition. Although majority of the calculations || ||,[|[L0| support the BKT 
phase transition curve starts from the point with the zero fugacity (or density) and the reduced temperature 0.25, we 
note that in the diagrams obtained from Monte Carlo simulations 0, the BKT phase transition curve starts from a 
slightly lower, but evidently different temperature. Here, this value of temperature, 0.25 corresponds to Coleman's 
critical point for the sine-Gordon (sG) continuum field theory jLl[ . 

The present disagreement on the phase diagram of the 2DNCCG is understandable. Actually, none of the above 
methods are exact. The renormalization group equations Q (1974) Q are valid just for low fugacity, so are Monte 
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Carlo simulations as reviewed in Ref. fL2| . Although an infinite subset of Feynman diagrams are re-summed up in 
Minnhagen's procedure [||, the generalized renormalization-group equations are also valid only for small fugacities 
when the temperature is lower than 0.25. Truly, the Gaussian approximation can work for both small and finite 
fugacities, but one does not know how good the results are. As for the generalized Debye-Huckcl-Bjerrum theory ||, 
it is of mean-field type intrinsically. 

The above analysis motivates us to propose a systematic, non-perturbative [], and, at the same time, convergent 
expansion scheme on the 2DNCCG to construct the phase diagram. This paper reports our effort of such an expansion 
which we call the non-perturbative sG expansion. Exploiting the equivalence 13 between the 2DNCCG and the sG 



field theory | |11|,[L4[ , we perform an expansion on partition function of the relevant sG field theory in Euclidean space 
to calculate the thermodynamic potential of the 2DNCCG up to the third order. The truncated result at the first 
order is equal to that obtained from the Gaussian approximation Jl0[ ] and, accordingly, the truncated results at higher 
orders are expected to provide systematic corrections to the Gaussian-approximation result. Thus, this scheme can 
afford an access to the correct phase diagram of the 2DNCCG. However, the higher-order corrections may give rise 
to additional divergences. Therefore, we treat, in the present paper, just the small-fugacity case for the second- and 
third-order corrections. In this case, taking a finite ultraviolet cutoff in the integration over the momentum gets rid 
of all the divergences appeared in the truncated expressions. Analysis of a characteristic parameter which is inversely 
proportional to the Debye screening length allows us to obtain phase diagrams of the 2DNCCG up to the third order. 
Although our study is limited to the small-fugacity effects of higher-order corrections, still, it provides systematic 
improvements towards the correct phase diagram. 

Based on the equivalence between the 2DNCCG and the sG field theory, the grand canonical partition function for 
a 2DNCCG consisted of point charges ±q Q |l(|(1997) becomes a functional integral over a real scalar field in the 
sG formalism, 

Z = Z- 1 J V[<j>]e~ I ' d2rn ^(r) ; (1) 
with Z = J V[4>] exp{J d 2 r\[V(j)(r)\ 2 }. The effective Hamiltonian density for the 2DNCCG, is given by 

n sG (r) = i[V0(r)] 2 - J cos[/30(f)], (2) 



where z is the electron-self-energy-renormalized fugacity, a an arbitrary length scale, and (3 = y with T the 

demensionless reduced temperature. Hereafter, we use <fi instead of <p{r). Calculating the partition function Eq.(l), 
one obtains the thermodynamic potential density Q, = — T\a[Z]/ J d 2 r, which is sometimes called free energy. 
Non-perturbative sG expansion, we intend to perform, starts with a modified Hamiltonian density, 

n sG (f,S)=H^ + S H D {r), (3) 

which is similar to what Refs. [p^l pTf did for the A0 4 field theory, where 5 is a parameter for the convenience of 
performing the expansion. Here, "H^Jf) = \<j){— V 2 + [i 2 )(f) and Hc(f) = —\p?$ 2 — |f cos[/3(/>]. When 6=1, the 
mass parameter fi cancels out, and, consequently, H s G(r, 5) is reduced to Eq.(2). Letting H S G(r, 5) take the place of 
Ti-scif) in Eq.(l), one can have the following modified partition function Zg, 

Z S = Z -% < e ^Jd 2 mo(f) ^ ( (4) 

with 

Z„ = J T>[</>] exp{ J d 2 r^(~V 2 + ^} (5) 



and 



< 



O > M = Z- 1 J V[fi]0 exp{ J d 2 r^{-V 2 + fi 2 )<j>} (6) 



Here, "non-perturbative" means that the coupling of the interacting system needn't be weak. 
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for an operator O. Expanding the exponential function in Eq.(4) into a Taylor series and, then, making a further 
expansion according to the series expression of the function ln(l + x), one has 



ln[Z 5 ] = - ln[Z ] + ln[Z M ] + ]T ^ ^-5 n < ( / d 2 rH D {^)) n >,] 1 

1=1 n=l ' J 

= - ln[Z ] + ln[Z M ] - 6 < J d 2 rU D (f) > M 

+S 2 \[< ( / d 2 rn D (r)) 2 > P -(< J d 2 rH D (r) > M ) 2 ] 

< ( / d 2 rH D (f) f >„ +~« | rf 2 rK fl (r) > M ) 3 
-i < J d 2 rH D {?) >,< ( J d 2 rH D {^)) 2 >„] + •••. (7) 

This is a power series in 6. Extrapolating the series to the case of 5 = 1 we can get the true thermodynamical potential 
density £1 which is independent of /i. However, it is impossible to calculate the above series exactly, and, therefore, 
one has to truncate the series at some order of 8 to obtain an approximation to CI. In the approximation, arbitrary a 
does not cancel out any more. Nevertheless, it can be fixed according to the "principle of minimal sensitivity" fL8|| , 
which means that, for a credible approximation, f2 should be insensitive to variations in \x. The general principle is to 
minimize the approximate result with respect to \x. Thus, we have a systematic tool to approximate Q order by order. 
Note that the above expansion is similar to the sG expansion in Ref. m there corresponds to fj, here) where A 
remains same for all orders in the expansion. However, we emphasize that, a crucial distinct point of our expansion is 
that the optimized fi changes from one order to the next. It is this point that ensures the present expansion convergent 
]T^ |. In fact, our expansion is just a generalization of the scheme proposed in Ref. JlTj]. We call it non-perturbative 
sG expansion. 

Now, we give a concrete calculation of fi up to the third order. In the right hand of Eq.(7), the first two terms arc 
easily calculated by employing the result of Gaussian functional integral |l9f] , and one can find the results in Ref. . 
As for the other terms, one can calculate them with the help of the external source technique of functional integration 
p9[ . In the calculation process, it is convenient to adopt the exponential form of cosine function, and the calculations 
can be performed as done in Ref. [[t0| (1990) and Ref. [^(J. Note that the averages of multi-cosine products can also 
be calculated in the same way in Ref. jl0| (1990) and Ref. Thus, a lengthy calculation leads to the following 

thermodynamical potential density fl, 



nr- 1 = 



f d2p 1 fix 1 / d2 P 1 r 2 _l 21 

-^ 2 / (1) [ M 2 ]-|exp{-i/3 2 / (1) [ M 2 ]} 

-^ (2 V] + i§/3V/ (2 V]ex P {-i/3 2 / (1) [M 2 ]} 



^(Jfexpi-^,)^]}^-^^^)^ 2 



•^ 6 / ( 3)[M 2 ] + ~^{-\i3 2 i {l) [AW 2 ^ki)[A - J/3V(/ (2 V]) 2 ] 

i(J) 2 exp{- / 3 2 / (1) [ M 2 ]}[/3 2 M 2 /( 2 )[^]^ 7 ^ j /3 4 ^^)[ A1 2 ] 



A- 



- (2k)\ 



/JV y 1 ; g2(2fc+l) ]J (2+2fc+l) ^2]] 



£(2* + 1)1 



4(!) 3 exp { -^ (1)[ , 2 ] }[ i ± ^r^ r^ -r^ ^ 

i,j=l y 1 ■> 

I \ " o2(2»+j+2fc-l) /-(2i-l ,j,2k) r 2i 
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i,J=l i,j,k=l 



where, 



2i-n 



,n-l ,2 n — 1 

/('"Vl = ( / dV)- 1 / d 2 r ld 2 r 2 d 2 r 3 ( \ i Pl * a , 2 )'( ' ^ P2 6 



and 



(2tt) 2 p 2 + /i 2 y 7 (2tt) 2 p 2 + /i 



JC-." V] = ( / rfV)- 1 / d 2 r ld 2 r 2 d 2 r 3 ( / ° g L , ■ 2 )' 



(2tt) 2 p 2 +/i 



(2tt) 2 p 2 +fi 2 ' 7 (2tt) 2 p 2 + ^ 2 ' • 

In the right hand side of Eq.(8), the first four terms represent the first-order result, next three terms the second-order 
correction and all other the third-order correction. 

Eq.(8) has many integrals and series summations. Obviously, many integrals in Eq.(8) are divergent. Introducing 
an ultraviolet cutoff in those integrals can remove the divergences. This is similar to what occurs in the Gaussian 
approximation []i"o| , ^i"|| . However, those series in Eq.(8) may still not be convergent for larger (3. In fact, this peculiarity 
is similar to the celebrated divergent problem of the sG field theory in the Minkowski space-time. For the (1+1)- 
dimensional sG field theory, Coleman's normal-ordering procedure H] can get rid of the divergences in the integrals 
in perturbation theory j2^| and Gaussian approximation |ll| , p3| ]. But, at any higher order of the perturbation theory 
for the range, An < 1 < 8ir, some series are still divergent. In the late 1970s and 1980s, several papers proposed a 
special treatment to the divergent problem of the sG field theory for the range An < (3 2 < 8n and performed further 
renormalization procedure to turn it finite p^j24| , in addition to Coleman's normal-ordering procedure. As afore 
mentioned, in this paper, we do not intend to treat the low-temperature divergent problem in Eq.(8). We note that in 
the Gaussian approximation, only linear terms in z survive, and the connection between the transition curves of the 
first-order and the BKT phase transitions involves only very small values of z. On the other hand, in Eq.(8), there 
exist terms with z 2 or z 3 , and only such terms with z 2 or z 3 adhere to series summations. Here, we retain only the 
linear terms in z, and neglect the higher-order terms of z in Eq.(8). This approximation allows us to compare the 
present result with the Gaussian approximation and observe the improvement directly. 

Retaining just the linear terms of z in Eq.(8) and carrying out the integrals with i as the ultraviolet cutoff, one 
can get the dimensionless reduced thermodynamical potential density £l = J7a 2 T _1 , 

Q = — ln[l + p 2 } - 2z(l + /2~ 2 )"^ 

87T 

1 p 2 Z(3 2 1 . __ 2 N-£l 

-(l+/i 2 ) K 



167T 1 + p 2 An 1 + p 

1 p2{1+2p2) < zp2 1 ra+^r^a + s, 2 )-^] ( 9) 



48tt (1 + p 2 ) 2 4tt (l + p 2 ) 2 " n ' L 2 V ' ' 16n' 

with p — [ia. In Eq.(9), the first two terms are just the Gaussian approximate result |[o]](PRE), the next two terms 
represent the second-order correction and the remaining terms the third-order correction. According to the "principle 
of minimal sensitivity" , p can be determined by minimizing Cl over p. That is to say, p should be satisfied with both 
the stationary condition ( Jp- = 0) and the stabilized condition ( (gjMa > 0). Note that /i is inversely proportional to 
the Debye screening length [js|J]~o|| , and, accordingly, p = implies that the gas becomes charge-binding and insulating, 
whereas p ^ signifies that the gas is conducting. Thus, we can obtain the phase structure of the 2DNCCG from 
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calculated values of p and Si. We have carried out a careful analysis of the results order by order. Fig.l is the T— z 
phase diagrams of the 2DNCCG of various orders constructed from the above calculation. 

In Fig.l, the dotted, dashed and solid curves correspond to the first-, the second- and the third-order results, 
repectively (As for the long-dashed line, see the explanation in the next paragraph). In the first order which is the 
Gaussian approximation, the T~z plane is divided into four regions : I, II, III and IV. The dotted curve between 
regions I and II is plotted from the stationary condition and the critical stabilized condition t^^y = 0. The dotted 
curve between II and III is obtained from Jl|p=o = il|p^o an d coincides with the solid curve in Fig.l of Ref. jl0| 
(PRE). The other two dotted curves which form the boundaries of region IV arc obtained through calculation of SI 
point by point and, hence, are rough, but qualitatively correct. In region I, p vanishes and the gas consists of dipolc 
pairs. In region II and III, every point have two values of p, p, = and p ^ 0, and, accordingly, the gas is in the 
insulator-conductor coexistence phase. But in II, the reduced thermodynamical potential f2 at p = is lower than Q 
at p ^ 0, and, in III, fl at p = is higher than Jl at p 7^ 0. In region IV, p and the gas is in the conducting phase. 
Analyzing the behaviour of the reduced thermodynamical potential f2, one can readily see that the dotted curve 
between regions I and IV corresponds to a continuous phase transition, i.e., the BKT phase transition, and on other 
dotted curves, first-order phase transitions occur. All of them meet at the tricritical point {T,z} = {0.25, l/167r}. 

An added feature in the second and third order results is the region V, which is separated from the region III by the 
long-dashed line |J Therefore, up to the second and third orders, the dashed and solid curves with the long-dashed 
line respectively divide the T—z plane into five regions: I, II, III, IV and V. Regions I, II, III and IV have the same 
physical meanings respectively with those for the case of the Gaussian approximation. As for region V, there exist 
two different non-vanishing p, and, accordingly, the gas consists in two different kinds of plasma phases and is a 
conducting coexistence phase. Note that the long-dashed line ends at the dashed curve for the second order and at 
the solid curve for the third order. 

In addition to the introduction of a new phase and the substantial extension of the insulator-conductor existence 
region, both the second- and the third-order corrections change the tricritical point significantly from that in the 
Gaussian-approximate phase diagram. The insulator-conductor coexistence region passes the end point of the BKT 
curve and end at the conducting phase. This peculiarity qualitatively agrees with the results from the Monte Carlo 
simulations on a triangular and a square lattice [Q and also from a generalized renormalization-group method 
We also note that the value of z at the lowest point of the dashed curve between regions II and III is 0.0238 and the 
corresponding z for the third-order case is 0.0255. These values show clearly a tendancy of approaching the value 
0.0291, which is obtained from the generalized renormalization-group method jfj]. 

Comparing our Fig.l with Fig.l in Ref. jl0]](PRE), one can see that the bottom of the dotted curve between regions 
I and IV, the BKT curve, in this paper is shifted to the right of the line T = 0.25, which is given in Ref. pO|](PRE). 
In fact, to get the the BKT critical line, Ref. pi (PRE) made a further approximation to Eq.(7) there. However, the 
present curve is obtained by directly checking the value of p from the Gaussian-approximation expression, without any 
further analytical approximation, and the second- and third-order corrections confirm this shift. We also note that the 
BKT critical curves obtained from the RG and generalized RG equations |^|-|6| start from the point {T, z} = {0.25, 0}. 
Nevertheless, We believe that the BKT critical curves in Fig.l do not conflict with the Coleman's phase transition 
point ^ and the RG results. Actually, only when the ultraviolet cutoff is enforced to infinity, the Coleman's phase 
transition occurs at (3 2 = 8ir |ll]]. The RG and the generalized RG equations take no such cutoff, whereas the 
present paper introduces a finite cutoff. When a finite cutoff is taken, the Coleman's phase transition does not occur 
and the sG continuum theory is well-defined at 1 — 87r and even over [ p5[ . 

This paper has performed a non-perturbative sG expansion on the 2DNCCG whose truncation at the first order 
is just the Gaussian approximation. We obtained the phase diagram of the 2DNCCG with the small-fugacity effect 
of the second- and the third-order corrections. Our result supports the prediction made by Minnhagen et al. || 
that the insulator-conductor coexistence phase ends at the conducting phase. Fig.l shows that the correction to the 
Gaussian approximation is important, and the smallness of the third-order correction presages fast convergency of 
the expansion in the present scheme. Therefore, we believe that our result for small fugacities is reliable within the 
framework adopted in the present paper, although, as reviewed in Ref. |h]](PRE), there is a source for uncertainty 
owing to the fact that the mapping between the 2DNCCG and the 2-dimensional sG field theory is exact only for the 
Euclidean sG field continuum without any finite cutoff . 



2 It is difficult to determine the left boundary of the region V. Here, the long-dashed line is given just to qualitatively indicate 
the existence of the region V. 

3 The Coleman's phase transition means that, for the 1 + 1-dimensional sG field continuum, its vacuum-state energy becomes 
unbounded from below when the coupling (3 2 increases to 8-7T, and hence the theory is ill-defined at that coupling. 
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In closing the paper, we intend to point out that it is straightforward to consider further higher order cor- 
restions. Moreover, Minnhagen's RG method || is based on a sG expansion, and therefore it is possible to generalize 
Minnhagen's method based on the expansion in the present paper. Perhaps, more interesting and meaningful, a 
complete treatment to Eq.(8), i.e., renormalization of the whole expression in Eq.(8) may yield a full phase diagram 
of the 2DNCCG. Finally, we intend to emphasize that the scheme here can be applied to other systems which is 
relevant to the sG model, for example, the backward-scattering model, the massive Luttinger model and Copper 
Benzoate p5|-p7|], and more general, even to other similar exponential-type interaction systems, such as, the massive 
sine-Gordon model p3| , 2D superconductor Q and finite-demensional site-disordered spin system |29| . 
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Fig. 1 Comparison among the phase diagrams of the 2DNCCG obtained in the Gaussian approximation, and with 
low-fugacity effects of second- and third-order corrections. The dotted, dashed and solid curves correspond to the 
Gaussian-approximation, second- and third-order results, repectively. The dotted curves divide the T—z plane into 
four regions: I,II,III and IV, and both the dashed and the solid curves with the long-dashed line yield an additional 
region V, besids regions I, II, III and IV. Region I corresponds to an insulating phase, regions IV and V to the 
conducting phases, and regions II and III to the insulator-conductor coexistence phases. Detailed explanation of this 
figure is given in the text. 
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